Power-law distributions for the areas of the basins of attraction on a potential energy 

landscape 
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Energy landscape approaches have become increasingly popular for analysing a wide variety of 
chemical physics phenomena. Basic to many of these applications has been the inherent structure 
mapping, which divides up the potential energy landscape into basins of attraction surrounding 
the minima. Here, we probe the nature of this division by introducing a method to compute the 
basin area distribution and applying it to some archetypal supercooled liquids. We find that this 
probability distribution is a power law over a large number of decades with the lower-energy minima 
having larger basins of attraction. Interestingly, the exponent for this power law is approximately the 
same as that for a high-dimensional Apollonian packing, providing further support for the suggestion 
that there is a strong analogy between the way the energy landscape is divided into basins, and the 
way that space is packed in self-similar, space-filling hypersphere packings, such as the Apollonian 
packing. These results suggest that the basins of attraction provide a fractal-like tiling of the energy 
landscape, and that a scale-free pattern of connections between the minima is a general property of 
energy landscapes. 
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The potential energy surface, which defines how the 
potential energy depends on the coordinates of all the 
atoms in a system, has a complex, multi-dimensional 
landscape In recent years, there have been inten- 
sive efforts to understand the behaviour of systems, such 
as proteins, supercooled liquids and clusters, in terms of 
features of these 'energy landscapes' Q. This research 
programme has led to important new insights into how 
proteins fold [H and the origins of the unusual dynamic 
properties of supercooled liquids [3, @ • 

Many of these applications rely on the 'inherent struc- 
ture' mapping introduced by Stillinger and Weber [(| that 
is illustrated in Fig. 1. It divides the energy landscape 
into basins of attraction surrounding the minima on the 
energy landscape, where a basin is defined as the set 
of points for which following the steepest-descent paths 
downhill from those points leads to the same minimum. 
The utility of the inherent structure approach is that 
it allows the behaviour of the complete landscape to be 
conceived in terms of the properties of the individual 
basins 0], which are themselves tractable to calculate. 
Hence, landscape-based descriptions of a system's ther- 
modynamics and dynamics can be obtained 

In many of the analyses of energy landscapes, the aim 
has been to understand differences in behaviour, e.g. pro- 
teins that are good or bad folders [1, Q , or supercooled 
liquids that show 'strong' or 'fragile' dynamics [ToL ITT1]. 
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in terms of differences in the energy landscape. However, 
surprisingly little is known about the more universal or- 
ganizing principles that are common to all such complex 
landscapes. For example, it is known that the number 
of minima scales exponentially with system size [l2j , and 
the distribution of minima as a function of energy is a 
Gaussian fl3| . but what is the nature of the division of 
the landscape into basins and the pattern of connectiv- 
ities between these basins? In particular, does the in- 
herent structure mapping lead to an equitable division 
of configuration space into basins, or an inequitable one 
where a small minority of the basins occupy a substantial 
majority of the space? 

Clues from recent work on the energy landscapes of 
small Lennard- Jones clusters [l4| [ll| perhaps suggest the 
latter scenario. The network of minima connected by 
transition states was found to be scale-free [Hj], that is 
the probability distribution for the degree (the number 
of connections to a node in the network) has a power-law 
tail. Such a feature has been found for many techno- 
logical, social and biochemical networks [l7], Gil, how- 
ever, unlike these other networks where the scale-free be- 
haviour arises from the dynamics of network growth [l6| , 
these 'inherent structure networks' are static. Hence the 
origin of the scale-free topology remains a puzzle. 

One potential answer is that the network structure re- 
flects an inequitable division of configuration space with 
the low-energy, high-degree minima having the largest 
basins of attraction [l9|, [20(] , because they are then more 
likely to have more transition states along the long basin 
boundaries. Indeed, a recently proposed model spa- 
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FIG. 1: (Color online) (a) A model two-dimensional energy 
landscape and (b) its associated contour plot illustrating the 
inherent structure mapping. In (b) the landscape has been di- 
vided into basins of attraction where the basin boundaries are 
represented by the thick lines, and the minima and transition 
states by points. 



tial, scale-free network based upon Apollonian packings 
2 ill has just such a correlation between area and degree 
22L l23l ]. A two-dimensional example of an Apollonian 
packing is illustrated in Fig. [2] Such a packing is gener- 
ated by the iterative addition of ever-smaller disks into 
the interstices of the packing until the space is filled. 
Consequently, the packing is fractal and self-similar [24[ • 
Importantly, the network of contacts between the disks 
is scale-free, and the properties of these 'Apollonian net- 
works' are very similar to the inherent structure networks 
[22I l23l ]. Thus, there seems to be a potential analogy 
between the hyperspheres in a high-dimensional Apollo- 
nian packing and the basins of attraction on an energy 
landscape. However, the starting point for this argument 
was the energy landscapes for clusters with less than 15 
atoms, and it is not clear whether the energy landscapes 
for these very small systems are representative. 

The fractal nature of the Apollonian packings is ev- 
ident from the distribution of the radii of the disks or 
hyperspheres making up the packing. This distribution 
is a power-law at small radii with the exponent related 
to the fractal dimension, df, where d — 1 < d f < d and 
d is the dimension of the space being packed [20 . This 
feature potentially allows the analogy between the Apol- 
lonian packings and the energy landscapes to be tested 
further. The distribution of the volumes of the hyper- 
spheres in an Apollonian packing scales as V~^ 1+df ' d \ 
thus leading to the prediction that, because of the high 
dimensionality of typical configuration spaces, the distri- 
bution of the hyperareas of the basins of attraction will 
follow a power-law with exponent —2, if the basins tile 
configuration space in an Apollonian-like manner. 

To test this prediction we need a method that can ob- 
tain this basin area distribution. To achieve this we study 
a transformed potential energy surface that is commonly 
used in global optimization and is particularly associated 
with the "basin-hopping" approach [2?], [H[ . This trans- 
formation involves assigning the energy at a particular 




FIG. 2: An Apollonian packing of a circle. Such a space- 
filling packing of disks is obtained iteratively, starting from 
an initial configuration of three mutually touching disks inside 
the circle. At each iteration disks are added inside each empty 
curvilinear triangle, such that each disk touches all three disks 
bounding the triangle. This process is repeated ad infinitum 
creating a space-filling fractal packing of the circle. Its fractal 
dimension is 1.3057 [251 ] . 

point in configuration space to that obtained after per- 
forming a local minimization from that point. It trans- 
forms the landscape into a set of steps, where each corre- 
sponds to a basin of attraction surrounding a minimum 
on the original potential energy surface. Thus, for the 
example landscape in Fig. [TJ the transformed landscape 
would consist of nine steps each having the energy of the 
corresponding minimum. Importantly, the probability of 
sampling a step during a simulation on the transformed 
landscape is proportional to its hyperarea [H, HH • 

Specifically, if we assume that the properties of the 
basins are uniquely characterized by the energy of their 
minima, the probability of being on a step with potential 
energy in the range E ± dE/2 in the canonical ensemble 
obeys 

^f n ns (£, T)dE cx n min (E)A(E) cxp (-E/kT) dE, (1) 

where f2 m i n (i£) is the number of minima with energy 
E ± dE/2 and A{E) is the average area of the basins 
of attraction surrounding minima with energy E. Hence, 
we can obtain both A(E) (to within a multiplicative con- 
stant) and the basin area distribution from p t ^ 1 f^ s (E,T) 
distributions obtained from simulations, if we first know 
Qmin(E). Methods to obtain Q, m i n (E) have already been 
developed [H, H3] , and involve inverting similar proba- 
bility distributions for the original potential energy land- 
scape. Namely, the probability of being in the basin of 
attraction of a minimum with potential energy in the 
range E ± dE /2 in the canonical ensemble is given by 

Pmin 

(E,T)dE cx tl m i n (E)Z v fo(E, T) exp (— E/kT) dE, 

(2) 
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where Z v ib(E,T) is the vibrational partition function of 
minima with energy E, which can be calculated, for ex- 
ample, using the harmonic approximation with frequen- 
cies obtained by diagonalization of the Hessian matrix of 
the relevant minima. 

Here, we apply this scheme [3ll | to three representative 
systems that have been much studied in the supercooled 
liquids community [32l . [33l l3~ij , namely a binary Lennard- 
Jones mixture with approximate composition A4B and 
Lennard- Jones parameters as given in Ref. l35l ; a model 
one-component glass-forming liquid interacting with the 
Dzugutov potential [36| ; and amorphous silicon modelled 
by the Stillinger- Weber potential [331 with a strengthened 
3-body term 34[. All three systems have 256-atoms and 
are modelled using periodic boundary conditions. 

For each system, we ran two sets of simulations. 
Firstly, constant-temperature molecular dynamics was 
performed and by regularly minimizing configurations 
generated along the trajectory we obtained p m - ln (E,T) 
and hence tl m m(E) (and a Gaussian fit to it). We chose 
temperatures well above the glass transition where equi- 
librium sampling of the liquid is easy to obtain. Secondly, 
Metropolis Monte Carlo simulations were performed on 
the transformed landscape, from which p t ^ 1 f^ s (E,T) and 
hence A(E) was obtained. As a check of the sampling, 
p t ^ 1 f^ s (E,T) distributions were collected at several tem- 
peratures. The resulting area distributions were in good 
agreement, confirming the reliability of the approach. 

In Fig. [Ha), A(E) is depicted for these three sys- 
tems. In agreement with our expectation that the deeper, 
more connected minima should have larger basins of at- 
traction, the basin areas decrease very rapidly with in- 
creasing energy in an approximately exponential man- 
ner. The probability distributions for the hyperareas of 
the basins of attraction are illustrated in Fig. E^b). It 
is apparent that the distributions follow an approximate 
power-law over the whole range of areas sampled, which 
is between 13 and 18 decades depending on the system. 
Closer inspection shows that the distributions begin to 
curve slightly downwards for larger areas, but that for 
the smaller basins the distributions very closely follow 
the prediction from the Apollonian analogy for 8 to 9 
decades for the binary Lennard- Jones and Dzugutov sys- 
tems, and 5 decades for amorphous silicon. 

These results have a number of important implications 
for the fundamental properties of energy landscapes, and 
the inherent structure mapping in particular. Firstly, 
they show that this mapping produces a very hetero- 
geneous division of configuration space with the deeper 
basins having much larger basins of attraction. Secondly, 
there appears to be a strong analogy between the way 
that basins of attraction tile the energy landscape, and 
hyperspheres fill space in an Apollonian packing. This 
similarity is exemplified by the power-law distribution of 
the basin areas at small A, which has exactly the ex- 
pected exponent over a very wide range of basin areas, 
and suggests that the basins divide configuration space 
in a fractal-like manner. Of course, the analogy must 
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FIG. 3: (Color online) (a) The dependence of the basin area 
on the energy of the minima, and (b) the cumulative basin 
area distributions for the binary Lennard- Jones (BLJ) and 
Dzugutov (Dz) liquids, and amorphous Si. To aid comparison 
between these systems the basin area and the energy are mea- 
sured with respect to the largest basin and the lowest-energy 
minima that have been sampled, respectively. Additionally, 
in (b) the power-law predicted by the analogy to Apollonian 
packings has been added for comparison [3||. 



break down at some sufficiently small length scale, be- 
cause the number of minima on an energy landscape, 
although large, is necessarily finite, whereas the Apollo- 
nian packings require an infinite number of hyperspheres 
to fill space. However, there is not yet any sign of this 
breakdown in the energy range that we have sampled. 
We should note that this energy range corresponds to 
the parts of the PEL sampled by a liquid, and in future 
work it would be interesting to try to use enhanced and 
biased sampling techniques to sample the basin area dis- 
tribution at both higher and lower energies. 

Thirdly, if, as for the Apollonian networks, the de- 
gree is a power-law function of the area [23[ (or equiva- 
lently the hyperlength of the basin boundary) the power- 
law distribution of basin areas is a signature that the 
underlying pattern of connections between the minima 
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is scale-free. Thus, our results suggest that this scale- 
free topology is not just specific to small Lennard-Jones 
clusters, but a more universal property of energy land- 
scapes. Indeed, there is evidence that the energy land- 
scapes of polypeptides also have a scale-free character 
[3!^ ]. Fourthly, the fractal- like character of the energy 
landscape provides a static explanation of the origin of 
the scale-free topology of inherent structure networks. 
However, this is not the end of the matter, since why the 
basins provide a fractal- like tiling of the energy landscape 



is still a puzzle, and one which we will explore in future 
work. 
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